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Abstract 

A method of exactly solving the master equation is presented in this letter. 
The explicit form of the solution is determined by the time evolution of a 
composite system including an auxiliary system and the open system in ques- 
tion. The effective Hamiltonian governing the time evolution of the composed 
system are derived from the master equation. Two examples, the dissipative 
two-level system and the damped harmonic oscillator, are presented to illus- 
trate the solving procedure. 
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The problem of open quantum systems has been around since the beginnings of quantum 
mechanics [1]. Important contributions to this general area have been made by researchers 
working in fields as diverse as cosmology [2], condensed matter [3], quantum optics [4,5], 
particle physics [6] , quantum measurement [7] , and quantum computation [8] . The problem 
can be described generally as interest in the effective dynamics of the open system surrounded 
by its environment. A formal framework to describe the effective dynamics of such an open 
system was set up in Ref.[9]. 

Generally speaking, interactions between a quantum system and its environment result 
in two kinds of irreversible effects: dissipation and decoherence. The first effect is due 
to the energy exchange between the system and its environment, whereas the second one 
comes from the system-environment interaction that does not change the system energy. A 
powerful tool to study the quantum dissipative system is the master equation, which can 
be obtained in Markovian limit [1,2,4]. This approximation is often very useful became it 
is valid for many physically relevant situations and became its numerical solutions can be 
found. As given by Gardiner, Walls and Millburn, Louisell in their textbook [4], the reduced 
density matrix p of the open system which is linearly coupled to its environment obeys the 
following master equation in Lindblad form [10] 
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Here, p{t) — p{t, K^, Gm) stands for the density operator of the system and p^nv denotes 
the density operator of the environment, X^ are eigenoperators of the system satisfying 
[Hq^X^ = ±hu!mX^, Hq stands for the free Hamiltonian of the system, and Am{Al^ 



are operators of the environment through which the system and its environment coupled 
together. Notice from eq.(l) that Gm should vanish at zero temperature T = 0, while 
Km should not if are indeed destruction operators of some kind. Some efforts such 
as the short-time expansion [11], the exact solution for some special kinds of the master 
equation[12], the small lose rate expansion[13] and the method of stochastic unravelings[14] 
have been made to solve (or solve it alternatively) the master equation. Unfortunately, this 
kind of master equation has not a general exact analytical solution yet. 

Here we shall construct an exact solution for the master equation to describe quantum 
dynamics of an open system with linear dissipation. The method is outlined as follows. After 
introducing an auxiliary system and mapping the density matrix of the original system to 
a pure state of the composite system, we obtain an effective Hamiltonian describing the 
evolution of the composite system from the master equation satisfied by the density matrix. 
Then the solution of the master equation can be obtained in terms of the evolution of the 
composite system by mapping the pure state back to the density matrix. As its applications, 
we present two examples to illustrate the solving procedure. The first example is a two-level 
atom coupling to a bosonic environment, whereas the second example consists of a single 
mode field in a cavity with linewidth k due to partial transmission through one mirror. 

We suppose that the Hilbert space of the system S is an A^-dimension Hilbert space 
spanned by, for instance, all the eigenstates of Hq. To begin with, we rewrite the master 
equation given in eq.(l) as the following general form 

i^^Hp-pH^ + tY: laLapLl (2) 

where H and La are well defined time-independent operators of the system, which may not 
be Hermitian generally and 7q, are real parameters determined by the environment. 

As the first step to solve this general master equation, we introduce an auxiliary system 
A which is the same with the original system we concern. That is to say, the original Hilbert 
space is extended to an A^^-dimensional Hilbert space, which is the Hilbert space of the 
original system S and the ancilla A. Let {|m5)|nA)} be an orthonormal and complete basis 
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for this composite system. 

A density matrix p of the system S*, whose matrix elements under the basis {l^ri^)} 
of the system are denoted as pmn = {TTT'\p\n)s, determines a pure bipartite state in the 
A^^-dimensional Hilbert space according to 

N 

iV'p) = E Pmn\ms)\nA). (3) 

m,n=l 

We note that this corresponding pure bipartite state is generally not normalized unless the 
initial state p of the system is a pure state. In fact every operator of the original system 
determines a pure state of the bipartite system in the same manner |T^. 

With the density matrix p evolving with time, the corresponding pure state \ipp) changes 
accordingly into another pure state 

N 

\Mi))= E Pmn{t)\ms)\nA). (4) 

m,n=l 

Because the evolution of the density matrix is governed by the general master equation (2), 
the evolution of the pure state li'pit)) is governed by the following Schodinger-like equation 

tdt\Mt))=H\ij,it))- (5) 
Here the effective Hamiltonian, which is generally not Hermitian, reads 

H = H - HA + i^laLaLAa, (6) 

a 

where operators H and of the original system are the same as given in Eq.(2) and 

operators Ha and are operators of the auxiliary system whose matrix elements are 
specified by 

{m\HA\n)A= {n\H^\m)s, (7) 

{m\LAa\n)A = {n\Ll^\m)s. (8) 

We note that the first two terms of the effective Hamiltonian H describe free evolutions of 
the system and the auxiliary system governed by H and —Ha and the third term describes 
an interaction between the system and the ancilla. 
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Because the effective Hamiltonian H of the composite system is time-independent, we 
can obtain formally the evolution of the pure state \ipp{t)) as 

\Mt))-e-'^'\^,). (9) 

^From the definition Eq.(4) of the pure state \ipp{t)), we obtain finally the time-dependence 
of the original density matrix as 

Pmn{t) = {ms\{nA\ipp{t)) = {ms\{nA\e-''''\^pp). (10) 

Since the effective Hamiltonian H changes its sign under the complex conjugate together 
with an interchange of the system and the ancilla, the Hermicity of the density matrix is 
preserved during the evolution. Moreover if the master equation (2) preserves the trace 
of the density matrix, solution given in Eq.(lO) has also trace 1 when the initial state is 
normalized. 

The generalizations of the above effective Hamiltonian method to the systems with in- 
finite many energy levels or with continuous spectra are straightforward. Also, the master 
equations with time-dependent ja and those master equations not posessing Lindblad form 
can be treated in the same manner, i.e., effective Hamiltonian can be obtained similarly. 
Thus the problem of solving the master equation becomes a problem of finding the time 
evolution of an effective Hamiltonian. Although the biorthogonal basis can be used to deal 
with general non-Hermitian Hamiltonian, for some special cases we can evaluate the evolu- 
tion operator directly for finite-level system as illustrated by the first example and factorize 
the evolution operator directly for quadratic systems as illustrated by the second example. 

In order to gain further insight into the content of the solution presented above, we 
shall consider at first the following master equation describing a two-level atom coupled to 
a bose-mode environment [4] 

P = P] + J{2cx-pcT+ - pa+a- - a+a.p} (11) 



with 
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where bl^{bm) stands for the creation (annihilation) operator of the m-th mode of the envi- 
ronment, Q is the Rabi frequency, and (J^((T±) denote the Pauh matrices. The basis of the 
two- level system is chosen so that (J+\gs) = \es) where l^f^) and le^) stand for the ground and 
excited states. This master equation is obtained under the condition that the environment 
is in its vacuum state. 

The auxiliary system we shall introduce is another two-level system whose Pauli matrices 
are denoted as Tz{t±) which satisfy t+Iqa) — \^a), where {ca) denotes the upper level of the 
ancilla A. By rewriting the master equation above into the general form of Eq.(2), we obtain 
immediately H = {flaz — ija+a-)/2 and L = a-. According to the definitions (7-8) we 
obtain Ha — {^t^ + i7T+T_)/2 and La — t^- The effective Hamiltonian governing the time 
evolution of the bipartite system is therefore 

H ^ H - HA + i-fa-T-. (12) 

As mentioned above, eq.(12) is the effective Hamiltonian corresponding to the master eq.(ll), 
which is derived by neglecting the temperature effect of the environment. If we take the 
temperature effects into account, the effective Hamiltonian should be 

Ht^ H - Ha + II (iVcT+T+ + (TV + 1)(7_T_ - N) 

^{f(..-r.)-n(iv4)} 

+i7 |tV(7+t+ + (TV + 1)(7_T_ - ^^^1 
= //o + ilJ (13) 

where TV = (exp(Q//cBT') — 1)~^ is the Bose distribution. Further we have [Hq,J] — 
and J\es)\9A) = J\9s)\^a) — which are sufficient to make an explicit calculation of the 
evolution operator 

/ TV + (TV + l)e-(2^+i)^* 7V(1 - e-(2^+i)^*) \ 

^ g-int-(iv+i/2)7t Q Q 

2TV-M gjm-(jv+i/2)7t Q 

^(TV+l)(l-e-(2^+^)^*) TV + 1 + TVe-(2^+^)^V 

(14) 
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where the bases of the composite system has been arranged as jeg, e^), jesjg'A), \gs-i^Aj-i 
and 15^5, qa). For simphcity, we consider zero temperature Hamiltonian H to study the time 
evolution of the density matrix p{t) with N — 0. 

For a general initial state p(0) = T,i,j=g,ePij\is){js\ where pij = {is\p{0)\js), the cor- 
responding initial state of the bipartite system is |V'p(0)) = J2i,j=g,ePij\'^s,jA)- With this 
initial condition the final state of the composed system at time t reads 

IV'pW) = Peee"^*|e5,eA) + {pgg + pee{l - e~'^^))\gs , Qa) 

+Pege-"^-''''\es,gA) + Pgee-^+''''\9s,eA). (15) 

That is 

Pee{t) = peee'"'\ pgg{t) = Pee(l " e"^*) + Pgg, 

p^g{t) = e-**-^Ve„ Pgeit) = e-^^+^^Vpe- (16) 

The results show that the off-diagonal elements of the density matrix are damped-oscillation 
function of time, while the diagonal element Pee decay exponentially. If Q > 7, there are 
several Rabi oscillations in p^g (or Pge) with the time evolution, otherwise (oi' Pge) decay 
directly. Especially we consider a state p — |e5)(e5| as the initial condition of the density 
operator p, At time t we have 

{es\p{t)\es) = e-^*, {9s\p{t)\9s) = 1 - e"^, 

{es\p{t)\gs) = {gs\p{t)\es)=0. (17) 

It is well known that the element {es\p{t)\es) of the density operator p{t) represents the 
population of the system in its upper level le^). The results (17) show that the population 
of the upper level le^) decay exponentially with the time evolution, this coincides with the 
results given in most textbooks [4]. 

The second example presented here is a single-mode field in a lossy cavity. The density 
operator for that mode obeys the following master equation in the Schrodinger picture[4], 

p = — i[a;/a^a, p] + —{2apo) — o)ap — po)a), (18) 
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where k is the hnewidth of the cavity mode with frequency cuf. In most textbooks, the 
solution of the master equation is given in terms of diagonal matrix elements (n|p|n) in 
a stationary state. Given an initial condition for the density operator, the evolution of p, 
however, is more useful than the stationary solution. In contrast with the solution in P- 
representation[4], in what follows, we present a solution of the master equation in a number 
state (Fock state) basis. 

Comparing the above master equation with the general form eq.(2), we see that 
H = {ojf—iK,/2)a) a and L = a. The auxiliary system is another single-mode field whose anni- 
hilation and creation operators are denoted as b and respectively. By definitions eqs.(7-8) 
we have Ha — {ojf -\- iK/2)b^b and La — b. Hence the effective Hamiltonian governing the 
time evolution of the bipartite system is 

H ^ (^^j: -—)a^a- (uf + ^-^)b^b + iKab, (19) 

The time evolution operator U{t) corresponding the effective Hamiltonian (19) reads 

JJ(^^^ ^ g-i(a;/-f )atatgi(a;y + f )bt6tgsta6^ ^20) 

where gt — 1 — e"***. We consider a general initial state p(0) = J2mn PmnlnT) {"^1, where 
Pmn = {'m\p{0)\n) is the element of the density operator at t = 0. Under the Fock state basis 
of the bipartite system, since \ipp{0)) — J2m,n Pmn\'m')\'n)b, with \n)fj satisfying b^b\n)fj — n\n)b 
we arrive at 

^ V(p + m)!(g + m)! ^ 

X 2^ Qt Pp+m,q+m- l^J-J 

m=o y'"'-) 

It is obvious that all elements of the density operator are damped-oscillation function of 
time, the decay rate K,{p + q) depends on the sum of p and q, which are initial condition- 
independent. Furthermore, we notice that the diagonal elements of the density operator, 
which represent the population in the corresponding state decay with different damped rate, 
the higher the state, the faster the decay. When the single-mode is initially at a thermal 
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state p(0) = e with the partition function as its normahzation, we find that after time 
t the system is still at a thermal state p{t) — e'^^*^"-^"- with 

(5{t) =P + Kt + ln(l - e-^(l - e"'^*)). (22) 

Based on equation (21), for the bipartite system in the initial state corresponding to p(0) 

we have 

oo 

Pit) = E AmmO)Alit) 

m =0 

I „m 

Amit) = W^e-(^'"/+'^/^)"^"*a"^, (23) 
V m! 

which shows explicitly that the evolution is completely positive and trace preserving. There- 
fore there exist an environment and an evolution U (t) of the system and the environment 
such that 

p{t) = TrE{U{t){p{0) ® |0)(OU)[/t(i)) (24) 

satisfies the master equation. It should be emphasized that the environment is exact. For 
the simple case of a single-mode field in lossy cavity, the environment is another single-mode 
field and the evolution is 

U{t) = e-^'"^*"^"e^*("^''-''^") (25) 

with cos^t = e"*^*/^. In fact Ajn{t) — {mb\U {t)\Ob) ■ This provides another Hamiltonian 
approach to this kind of problem: two-mode field interacting with one another with a time- 
dependent interaction. 

In summary, we propose a method to approach solving the master equation exactly. For 
this end, we first of all introduce an auxiliary system, which has the same dimension as 
the system that we are interested in. The original system and the ancilla interact on each 
other, and the Hamiltonian which governs the time evolution of the composite system are 
determined by the master equation. In this sense, the solution of the master equation might 
be computed through the Schrodinger equation of the composed system. We would like to 

9 



note that whether we can obtain the solution of the master equation in an explicit form 
or not depends on the form of the effective Hamiltonian, but the method presented here 
provides a new approach to the exact solution of the master equation. 

This work was supported by Chinese postdoctoral Fund via Institute of Theoretical 
Physics, Academia Sinica, and by K. C. Wong Education Foundation, Hong Kong. 
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